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Dynamical localization, i.e. the absence of secular spreading of a quantum or classical wave packet, is usually 
associated to Hamiltonians with purely point spectrum, i.e. with a normalizable and complete set of eigenstates, 
which show quasi-periodic dynamics (recurrence). Here we show rather counter-intuitively that dynamical 
localization can be observed in Hamiltonians with absolutely continuous spectrum, where recurrence effects 
are forbidden. An optical realization of such an Hamiltonian is proposed based on beam propagation in a 
self-imaging optical resonator with a phase grating. Localization without recurrence in this system is explained 
in terms of pseudo-Bloch optical oscillations. © 2015 Optical Society of America 
OCIS codes: 000.1600,140.4780, 050.2770 


Introduction. Diffusion, localization and recurrence 
phenomena play a crucial role in several fields of sci¬ 
ence [1]. Optics has provided on several occasions an 
ideal framework to investigate transport and localiza¬ 
tion effects earlier predicted in quantum physics [2,3]. 
In systems described by a Schrodinger-like wave equa¬ 
tion, general relations between the temporal dynam¬ 
ics and the spectrum of the underlying Hamiltonian 
have been disclosed (see [4-10] and references therein). 
Roughly speaking, the absolutely continuous spectrum 
case corresponds to transport and diffusive phenomena, 
the pure point case typically corresponds to the absence 
of transport (localization) and recurrence phenomena, 
while the singular continuous case corresponds to inter¬ 
mediate transport behavior. Hamiltonians with a purely 
point spectrum always show almost periodic dynamics 
owing to the quantum recurrence theorem [11], whereas 
in systems with purely absolutely continuous spectrum 
recurrence effects are always prevented as a consequence 
of the Riemann-Lebesgue theorem [4,5]. In case of pure 
point spectrum, in addition to recurrence one generally 
observes dynamical localization^ i.e. an initially localized 
wave packet does not secularly spread, whereas in the 
absolutely continuous case dynamical delocalization is 
commonplace. Prominent examples of localization and 
recurrence associated to a point spectrum are the phe¬ 
nomena of Anderson localization and its ramifications, 
and of Bloch oscillations, which have been extensively in¬ 
vestigated in optics [12-29]. While the above rule is very 
common, it is not universal one. For example, dynamical 
delocalization in systems with pure point spectrum can 
be found in certain random Hamiltonian models [30]. 

In this Letter we present an example of a quan¬ 
tum Hamiltonian H with absolutely continuous spec¬ 
trum that shows dynamical localization, and suggest 
an optical implementation of this model. The Hamilto¬ 
nian H describes a continuous extension of the effective 
Wannier-Stark Hamiltonian in a single band of a crys¬ 
tal [31], i.e. it is of the form H = T{px) + Fx where 
the kinetic term T a periodic function of the momentum 


operator px = —idx with period 27r/d whereas x is a 
continuous spatial variable. While in the discrete limit, 
i.e. by taking x = nd with n = 0, ±1, ±2,... the Hamilto¬ 
nian H has a pure point spectrum (Wannier-Stark lad¬ 
der), showing localization and quantum recurrence (the 
famous Bloch oscillations), in the continuous case peri¬ 
odicity is lost according to the Riemann-Lebesgue the¬ 
orem, however dynamical localization is preserved. We 
show that such a dynamical behavior, that we may re¬ 
fer to as pseudo-Bloch oscillations, can be realized in a 
self-imaging optical ring resonator with an intracavity 
phase grating. A transition from localization to delocal¬ 
ization by breaking the exact self-imaging condition of 
the resonator is also discussed. 

Quantum Hamiltonian and pseudo-Bloch oscillations. 
Let us consider a quantum system described by the one¬ 
dimensional Schrodinger equation idti^{x,t) = Hfj{x,t) 
with a Stark Hamiltonian given by (ft = 1) 

H = T{px) + Fx, (1) 


where Px = —idx is the momentum operator, F is the 
external force and T{px) is the effective kinetic energy 
operator. When the kinetic energy operator is the one of 
a free particle, i.e. for T{px) = the energy spectrum E 
of H is purely absolutely continuous —oo<E<oo with 
improper eigenfunctions \(f)E{x)) given by shifted Airy 
functions. According to the Riemann-Lebesgue theorem, 
the return probability, defined as 


Pr{t) = \{^p{x,0)\^p{x,t))f 


/ 


dx2p*{x, 0)2p{x, t) 


2 


( 2 ) 

decays to zero. Moreover, dynamical delocalization oc¬ 
curs. In fact, by changing the reference frame from rest 
to an accelerated one, the Stark Hamiltonian p^ + Fx 
is equivalent to the one of the free particle without the 
force (see, for example, [2]). Here we consider the case 
where T{q) is a period function of q with period 27r(d, 


A9) = E Ti expQldq) (3) 
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with T-i = T^*. Hamiltonians with a periodic kinetic 
energy operator are generally introduced in solid-state 
physics as an effective Hamiltonian to describe single¬ 
band electron wave packet dynamics in slowly-varying 
external fields [31]. In this case T{q) is the band dis¬ 
persion curve, and the discretization x = nd {n = 
0, ±1, dz2, dz3, ...) at localized Wannier sites must be ac¬ 
complished [31]. As is well-known, with such a discretiza¬ 
tion H has a pure point spectrum with equally-spaced 
energy levels (Wannier-Stark ladder) and energy separa¬ 
tion ujb = Fd] in real space an initially localized wave 
packet undergoes a periodic motion, so-called Bloch os¬ 
cillations, and the return probability Pr{t) is a peri¬ 
odic function with the Bloch oscillation period Tb = 
271 jujB = 271/{Fd). Here we consider the case where the 
spatial variable x is continuous rather than discretized. 
In such a case, it can be shown that H has a pure ab¬ 
solutely continuous spectrum —oc<E<oc with im¬ 
proper (non-normalizable) eigenfunctions |0£;(x)) given 
by 

4>e{x) = Y,-nd- E/F) (4) 

n 

where we have set 

( 5 ) 

Note that \(I)e{x)) are not normalizable and satisfy the 
usual relations of improper eigenfunctions {(I)e'\4^e) = 
6{E — E'). The most general wave packet solution to the 
Schrodinger equation idt^ = H^p is given by 

V^(x,t) = / dEju{E)(l)E{x) (6) 

where the spectral measure ju{E) is determined by the 
initial wave packet distribution '0(x,O) according to 
/i(^) = J (ix(/)|;(x)'0(x, 0). The return probability Pr{t) 
then reads Pr{t) = |/ dE\fi{E)\‘^ ex.p{iEt)\‘^ ^ which de¬ 
cays to zero as t ^ oo according to the Riemann- 
Lebesgue theorem: hence when the spatial variable x is 
continuous rather than discrete, the dynamics is not pe¬ 
riodic. Using Eq.(4), Eq.(6) can be cast in the form 

'0(x, t) = G(x, t) ex.p{—iExt) (7) 

where we have set 

G(x, t) = F 0i0p^'ip{x — ld-\- nd^ 0) exp(iF(i/t) (8) 

n,l 

Note that G(x, t) is a periodic function of t with period 
Tb and G(x,0) = '0(x,O). Since \ip{x^t)\'^ = |G(x,t)p 
is periodic, dynamical localization is obtained. However, 
contrary to ordinary Bloch oscillations the dynamics is 
not periodic [this is due to the factor ex.p{—iExt) in 
Eq.(7)] and the return probability decays to zero. We 
can refer to such a dynamical behavior, i.e. dynamical lo¬ 
calization but absence of recurrence, to as pseudo-Bloch 
oscillations. An example of pseudo-Bloch oscillations in 


an optical resonator will be presented in the following. 
Einally, let us extend the analysis by considering the 
Hamiltonian 

H = T(pa;) + epl + Fx, (9) 

where T(g) is periodic and defined by Eq.(3) and an addi¬ 
tional non-periodic kinetic term of strength e is included. 
Eor e 7 ^ 0, it can be shown that pseudo-Bloch oscilla¬ 
tions are smeared out and dynamical wave packet delo¬ 
calization arises. In fact, the solution to the Schrodinger 
equation with Hamiltonian H defined by Eq.(9), from 
time tn = uTb to time tn-\-i = (n + I)Te, is formally 
given by 'ip{x, t^+i) = exp(—iTe77)'0(x, t^). By repeated 
use of the Baker-Campbell-Hausdorff formula and af¬ 
ter some cumbersome calculations, it can be shown that 
' 0 (x,tn+i) = exp(—iTe^e//)' 0(^5 ^n), where we have set 

Heff = epl + Fx (10) 

Such a result, which is exact regardless of the strength of 
e, shows that the quantum evolution of the Hamiltonian 
(9), mapped at times tn = uTb multiples than the Bloch 
oscillation cycle, is equivalent to that described by the 
effective Hamiltonian (10). Since Heff is equivalent to 
the free-particle Hamiltonian after a change of reference 
frame, it shows dynamical delocalization. 

Optical realization. Transverse light dynamics in op¬ 
tical resonators can provide a fertile testbed to emulate 
continuous quantum mechanical Hamiltonians with en¬ 
gineered kinetic energy operator [32-34]. An optical im¬ 
plementation of the continuous quantum Hamiltonians 
(I) or (9) with periodic T{px) can be obtained by con¬ 
sidering the optical resonator shown in Eig.I. It consists 
of a (near) self-imaging ring cavity of length L in so- 
called 4/ configuration, comprising four equal converg¬ 
ing lenses of focal length /, one phase grating (PG) and 
one prismatic phase mask (PM) in the geometrical set¬ 
ting of Eig.I. The total cavity length is L = 8/ -h /, with 
I f and / = 0 for exact self-imaging. The transmission 
functions of the PG and of the prismatic PM are given 
by ti{x) = ex.p[—iipi{x)] and t 2 = ex.-p[—iip 2 {x)], where 
ipi{x a) = (•pi{x) is the periodic phase change intro¬ 
duced by the PG with spatial period a and y72{x) = iakx 
is the phase change introduced by the prismatic PM 
(/c = 27r/A is the optical wave number and a is the angle 
of the prismatic surface). Light propagation inside the 
optical ring can be readily obtained by application of the 
generalized Huygens-Eresnel integral (see e.g. [32-34]). 
Assuming one transverse spatial dimension x and indi¬ 
cating by p7m{x) the envelope of the intracavity field at 
the reference plane 7 in the cavity at the m — th round 
trip, the following recurrence relation can be derived 

'Pm+l{x) = i>it>2'D3'prn{x) ( 11 ) 

where we have set Vi = exp(—X >2 = 

exp[—i(^ 2 (^)], and the operator P 3 is defined via the rela¬ 
tion 7 ) 3 /(x) = (Tnf{x-\-nd)^ with d = Xf /a and with 


2 




lens/ plane Y lens/ 


Fig. 1. (Color online) Schematic of the optical ring res¬ 
onator that realizes pseudo-Bloch optical oscillations. The 
cavity operates in the exact self-imaging condition for / = 0 . 
The resonator comprises four lenses (focal length /), a phase 
grating and a prismatic phase mask, with transmittances 
ti{x) = exp[—i(^i(x)] and t 2 {x) = eyip[—i(p 2 {x)] respectively. 
The total cavity length is L = 8 / + /. Planes 7 and 7 ^ are 
Fourier conjugate planes. 

coefficients given by the Fourier series ex.p[—i(pi{x)] = 
(Jn exp(27rmx/a). In writing Eq.(ll), we neglected 
cavity losses and assumed that a single axial mode of 
the cavity is excited [32,34]. Operator Vi on the right 
hand side of Eq.(ll) describes beam diffraction for an 
effective length / (/ = 0 and Vi = X for an exact self¬ 
imaging cavity), V 2 accounts for the additional phase 
acquired after transmission through the prismatic PM, 
whereas P 3 describes the effect of the PG placed in the 
Eourier plane (see Eig.l). To clarify the analogy be¬ 
tween Eq.(ll) and the Hamiltonian (9), let us consider 
the limits |(/i, 2 | ^ tt and / ^ 0 , so that the operators 
^ 1 , 2,3 slightly deviate from the identity operator X. In 
this case, using standard methods (see e.g. [32-34]) af¬ 
ter first-order expansion and continuation of the round 
trip number (m ^ t) from Eq.(ll) one can derive the 
following master equation that describes transverse light 
evolution in the resonator 

n 

where t is the temporal variable in units of the cav¬ 
ity round trip time, F = ka, e = —//(2/c), and are 
the Eourier coefficients of the phase (/i(x) [35]. Note 
that Eq.( 12 ) is precisely the Schrodinger equation with 
Hamiltonian (9). Hence for exact self-imaging condition 
(/ = 0 ) the beam evolution over successive round trips 
in the ring resonator undergoes pseudo-Bloch oscilla¬ 
tions and shows dynamical localization, whereas smear¬ 
ing out of pesudo-Bloch oscillations and delocalization 
occurs when the resonator length is detuned from exact 
self-imaging. An example of localization without recur¬ 
rence arising from pseudo-Bloch oscillations under ex¬ 
act self imaging condition is shown in Eig.2. The fig¬ 
ure shows the evolution of an initial Gaussian beam dis¬ 
tribution in the ring at successive round trips, as ob- 



round trip number t round trip number t round trip number t 


Eig. 2. (Color online) Pseudo-Bloch oscillations in the ring 
resonator under exact self-imaging condition (/ = 0 ). (a) and 
(b) show in pseudo color maps the numerically-computed evo¬ 
lution of the intensity [panel (a)] and real part [panel (b)] 
of the intracavity held '/(x,t) at plane 7 at successive round 
trips t in the ring. Parameter values are given in the text. The 
initial condition is the Gaussian held '0(x,O) oc exp(— 
with spot size re = 150 /xm. (c) Evolution of the correlation 
function Pr(t). 


(a) (b) 
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Eig. 3. (Color online) Smearing out of pseudo-Bloch 
oscillations in the ring cavity under a slight deviation 
from self-imaging condition {1 — 2 mm). Other parame¬ 
ter values are as in Fig.2. (a) Numerically-computed evolu¬ 
tion of the intensity of the intracavity held ^{x,t) at plane 
7 , and (b) of the beam spot size, dehned as [J dx{x — 
(x))^|'0(x, t)|^/ f dx|'0(x, t)|^]^/^. The secular growth (on av¬ 
erage) of the spot size at successive Bloch oscillations cycles 
is the signature of dynamic delocalization. 

tained from the propagation map ( 11 ), for a sinusoidal 
PG (pi{x) = 6 cos(27rx/a) and for parameter values / = 3 
cm, A = 633 nm, 6 = 0 . 2 , a = 180 /xm, a = 0.2 mrad 
and L = 8 / = 24 cm. Figures 2(a) and (b) show the 
evolution of the intracavity held intensity |' 0 (x,t)p and 
of the real part of the held Re[' 0 (x,t)], respectively, at 
plane 7 ; in (c) the evolution of the correlation function 
Pr{t) = 1/ dx/^*(x, 0 )/^(x, t)p (the analogue of the re¬ 
turn probability) is also depicted. Note that, while the in¬ 
tensity distribution undergoes a periodic dynamics with 
period Tp = 2ii/{Fd) = a/{af) 30 round trips, the 
dynamics is not periodic and the held does not repro¬ 
duces itself, as one can see from Figs.2(b) and (c). 
Smearing out of pseudo-Bloch oscillations and dynamical 
delocalization are observed when the ring deviates from 
the self-imaging condition. As an example. Fig.3 shows 
the beam dynamics for the same conditions of Fig.2, ex¬ 
cept for an ehFective diffraction length I = 2 mm in the 
ring. In this case one clearly observes that the oscilla- 
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tory dynamics of Bloch oscillations is superimposed to a 
transverse acceleration of the optical beam [Fig.3(a)] and 
a spreading of the beam size on average [Fig.3(b)]. Such 
a behavior arises from the fact that the beam dynam¬ 
ics mapped at successive Bloch oscillation cycles is de¬ 
scribed by the effective Hamiltonian (10). Given the form 
of Hef f , effective diffraction and transverse index gradi¬ 
ent are responsible for spreading and transverse beam 
acceleration. 

To experimentally observe pseudo-Bloch oscillations, 
one can consider the passive ring resonator of Fig.l 
driven by an holding beam coupled through one of the 
four mirrors of the cavity (see also Ref. [34]). Assuming 
that the frequency of the holding beam is in resonance 
with one of the cavity axial modes, the map ( 11 ) is re¬ 
placed by the following one 

1 pm+l{x) = 'DiT> 2 t>si;rn{x) + VTEm{x)-'^i)m{x), (13) 

where T ^ 1 is the transmittance of the coupling mirror 
and Em{x) is the spatial profile of the injected beam at 
the m-th round trip and at plane 7 {Em{x) is a slowly- 
varying function of m). Similarly, taking into account 
cavity losses and external beam injection the master 
equation ( 12 ) is modified as follows 

^ -e^^+Fx' 0 +^Tn' 0 (x+nd)-i^' 0 -hi\/TF;(x, t). 

n 

(14) 

Pesudo-Bloch oscillations can be observed when a 
Gaussian holding beam is switched off and the transient 
decay dynamics of the passive ring resonator (i.e. the 
transverse beam distribution over several round-trips 
during the transient decay) are monitored outside 
the cavity by a fast GGD camera, as discussed in 
details in Ref. [34]. As an example. Figs.4(a) and 
(b) show the numerically-computed evolution of the 
normalized intensity distributions of the field 'ip{x,t) at 
plane 7 after switching off the Gaussian holding beam 
E{x^t) = /(t) exp(—x^/rc^) with a decaying amplitude 
fit) = ( 1 / 2){1 — tanh[(t — to)/At]} [dashed curve in 
Fig.4(c)]; parameter values used in the simulations are 
A = 633 nm, / = 3 cm, a = 0.1 mrad, a sinusoidal PG 
^lix) = (5cos(27rx/a) with S = 0.3 and a = 180 jum, 
output coupling T = 1%, spot size of the holding 
Gaussian beam w = 180 /im, and switch off time 
At = 8 (in units of round-trip time). In Fig.4(a) exact 
self-imaging of the cavity (L = 8 / = 24 cm) is assumed, 
whereas in (b) a slight deviation from self-imaging by 
I = 1 mm is assumed. The solid curve in Fig.4(c) shows 
the evolution of optical power available outside the 
resonator during the switch-off dynamics, normalized 
to its continuous-wave value before switching off the 
holding beam. Note that, before the holding beam is 
switched off, a stationary intracavity field distribution 
is established, which deviates from the Gaussian one 
because of the effects of PG, transverse index gradient 
and [in (b)] diffraction. As the holding beam is switched 


off at round trip time to = 200 , the dynamical behaviors 
of beam intensity in Figs.4(a) and (b) become similar 
to those shown in Figs.2(a) and 3(a), respectively. Note 
also that, since the Bloch oscillation period is smaller 
than the photon lifetime in the cavity, a few Bloch 
oscillation cycles can be detected during the decay 
dynamics. 



round trip number t round trip number t round trip number t 


Fig. 4. (Color online) Pseudo-Bloch oscillations (color maps 
of normalized intensity distribution at plane 7 ) in the ring 
resonator with a Gaussian holding beam for (a) exact self¬ 
imaging, and (b) for a detuned cavity (/ = 1 mm). Parameter 
values are given in the text. Panel (c) shows the behavior of 
normalized optical power (solid line) and amplitude /(t) of 
the Gaussian holding beam (dashed line) during the switch- 
off time. 

Conclusions. In conclusion, the spreading properties 
of quantum or classical wave packets described by a 
Schrodinger-like equation are known to be closely re¬ 
lated to the spectral properties of the Hamiltonian. An 
absolutely continuous spectrum is generally associated 
to dynamical delocalization and absence of recurrence, 
whereas a pure point spectrum generally corresponds to 
localization and recurrence. However, such rules are not 
universal. Here we have shown rather counter-intuitively 
that dynamical localization can be observed in Stark 
Hamiltonians with purely continuous spectrum, and in¬ 
troduced the novel phenomenon of pseudo-Bloch oscil¬ 
lations. An optical realization of the quantum Hamil¬ 
tonian, based on transverse beam dynamics in a ring 
resonator, has been suggested. Our analysis indicates 
that light dynamics in optical resonators can provide a 
testbed for the observation of other unusual localization 
and diffusion processes. For example, by proper engi¬ 
neering the power law of effective diffraction in the res¬ 
onator [33], regimes like quasi-absence of diffusion pre¬ 
dicted in certain quantum models [ 8 ] could be observed 
in optics. 
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